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Abstract. We discuss interfaces in spin glasses. We present new theoretical results 
and a numerical method to characterize overlap interfaces and the stability of the spin- 
glass phase in extended disordered systems. We use this definition to characterize the 
low temperature phase of hierarchical spin-glass models. We use the Replica Symmetry 
Breaking theory to evaluate the cost for an overlap interface, which in these models is 
particularly simple. A comparison of our results from numerical simulations with the 
theoretical predictions shows good agreement. 
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1. Introduction 

The study of interfaces is often used to investigate the stabihty of ordered coexisting 
phases. Low temperature ordered phases are spatially homogeneous, however, interfaces 
can be induced by forcing heterogeneous boundary conditions [1]. For example, in the 
case of phases characterized by different values of an order parameter, one can fix the 
order parameter at the boundaries along one selected direction to two different values. 
For stable phases the induced interface costs an amount of free-energy proportional to 
a power of the interface size. At the lower critical dimension this power vanishes and 
long-range order becomes unstable. 

In spin glasses, due to the lack of a physical order parameter allowing to distinguish 
the different phases, this definition of an interface does not seem of practical use. 
Different proposals have been put forward to define interfaces in spin glasses. A 
common procedure [21 El H] , involves the comparison of free-energies between systems 
with periodic and antiperiodic boundary conditions, respectively. This leads to a 
definition of an interface exponent, from the relation: |Fper — -Fanti-per| ~ which 
plays a fundamental role in the droplet picture of spin glasses [U [5] . Measures of have 
been used to perform a high-precision estimate of the spin-glass lower critical dimension 
Dc ~ 2.5 in Ref. [6] (defined as the value of D for which the interface exponent 9 
vanishes). An alternative definition [7] uses the free-energy cost that is needed to impose 
a spatial heterogeneity in the order parameter, which in spin glasses is the overlap, i.e. 
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a measure of the similarity between two equilibrated configurations of the same system 
[8]. One needs to consider two identical copies of the system constrained in such a 
way to have their mutual overlaps on the boundaries fixed to some preassigned values. 
This definition is well suited to test an important property emerging in the Replica 
Symmetry Breaking (RSB) of the spin-glass phase: while the probability distribution 
p{q) for the overlap is broacjjl, reflecting the presence of many different and almost 
degenerate equilibrium phases, one has that for any couple of equilibrated configurations 
their mutual overlap is spatially homogeneous on a large scale. Two configurations with 
a given value of the overlap in a given macroscopic portion of space will have the same 
overlap everywhere. 

Starting from this observation a pair of identical systems subject to "twisted 
overlap" boundary conditions was considered in |7|]. Two different values qi and q2 
of the overlap were imposed on the two boundaries along a selected direction, chosen 
among the ones with non-zero probability density in the p{q). It was then argued, that if 
RSB is present, the boundaries induce a smooth overlap profile in space that interpolates 
between the two values. Above the lower critical dimension Dc for spin-glass order this 
has a diverging free-energy cost, and one can expect the following scaling with L and 

111 - <l2\- 



where the exponent 9q = D — Dc is positive for D > Dc. A mean-field calculation of 
bg and the interface exponent 9g gave the non-trivial values of 6^ = 5/2 and Dq = 5/2. 
Remarkably, while measuring a quite different property, the resulting value of the lower 
critical dimension is in very good agreement with the one estimated by the periodic- 
antiperiodic boundary condition method [HE]. 

Unfortunately, despite its potentially informative content on the nature of the spin- 
glass phase, the definition of |7] was never used in numerical studies of the spin-glass 
phase. This is probably due to the difficulty of imposing values of the overlap at the 
boundaries. 

In this paper we start from the general observation, in fact not at all specific to spin 
glasses, that it is possible to induce interfaces with an alternative procedure. One can 
divide the physical system into two contiguous halves and impose fixed values to the bulk 
order parameters in each of the two halves. This procedure, theoretically equivalent to 
fixing heterogeneous boundary conditions, is better suited for numerical investigations, 
since the resulting free-energy cost can be related to the probability of a spontaneous 
fluctuation in an unconstrained system. 

In the case of spin glasses we have to consider two real replicas and impose values 
qi and q2 of their mutual overlaps in two contiguous half spaces. This way of imposing 
heterogeneity in the system induces an overlap profile in space equal to the one induced 
by fixing the overlap on the boundaries to suitable values q[ and Oiig then has a 

I This is in contrast to the droplet picture of the spin-glass phase where the distribution p{q) - although 
having noticeable tails for finite systems - reduces to two delta functions in the thermodynamic limit. 
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free-energy cost AF(qi, 52) of the form ([T]) with the same exponents 9q and The free- 
energy cost for imposing the different overlaps is related by Boltzmann's relation to the 
probability that a spontaneous fluctuation of a couple of unconstrained systems produce 
the values qi and q2 of the overlaps in the two half spaces: PL{qi,q2) = e-^^^(<?i''?2), which 
is a large deviation formula. This quantity has the advantage that it is easily accessible 
in numerical simulations and therefore we base our analysis on it. 

In this paper we focus our attention on overlap interfaces in Dyson-like hierarchical 
spin-glass models. 

Hierarchical models without disorder [9] have played an important role in the 
theoretical and mathematical understanding of critical points [ini IH]. In the 
ferromagnetic Dyson model one can write exact renormalization group equations that 
involve the iteration of a function of a single variable (for a review see [l2])- Hierarchical 
spin-glass models of the Dyson kind have not to our knowledge been considered in the 
literature and are in our opinion very attractive as they could allow for an analytical 
study of non-mean-field disordered critical points and nontrivial low temperature phases. 
These models provide a hierarchical version of spin glasses with power-law interactions 
introduced in [13]. One-dimensional models with power-law interactions have recently 
received a great deal of attention as test grounds for theoretical ideas about finite-range 
spin-glass phases [m [151 [HI [13 [H]- Here, following Dyson, we propose to use the 
tree topology as a further simplification. We choose to start the study of these models 
focusing on interface properties, for which we can get particularly simple theoretical 
predictions in the RSB framework. We study the system through the replica method, 
deriving a recursion equation relating the replica partition functions at the different 
levels of the hierarchy. This relation, which in principle codes for all the thermodynamic 
properties of the system is analyzed in a self-consistent way to obtain the interface free- 
energy. We test our prediction in numerical simulations, finding good agreement. 

The organization of the paper is the following: in section [2] we discuss the definition 
of the interfaces and how in principle interfaces can be evaluated and some expectations 
based on RSB theory. Section [3] is devoted to the definition of the models we use. In 
section H] we discuss possible scenarios and expectations and in section [5] we sketch the 
theoretical derivation of the interface free-energy cost in hierarchical models. Section 
M discusses the numerical simulations. Finally, section [7] concludes the paper. In the 
appendix we discuss the numerical characterization of the critical point. 

2. Generalities 

In this section we define the probability distribution of overlaps in contiguous half spaces 
and discuss how it can be computed in principle in the large deviation regime. For 
definiteness we consider a spin-glass system with Ising spins Si = ±1 on a set of indices 

§ This is strictly true for a couple of systems with open boundary conditions where the first procedure is 
naturally defined. The second procedure has the advantage of being defined also for periodic boundary 
conditions. 
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2 G A with elements, that we divide in two "half spaces" Ai and A2 with N/2 elements 
each. For two spin configurations S and S' one can define the partial overlaps Qr{S, S') 
in each of the half spaces A^, r = 1,2 as Qr{S,S') = i/{N/2)J2i,^A^ SiS[. The joint 
probability distribution function (PDF) of the two overlaps, for fixed value of the system 
size L and quenched disorder J, is given by 

S,S' 

where Z denotes the partition function. From this relation one gets the usual PDF of 
the total overlap by simple integration. One can consider the average over the disorder 
of ([2]), however, here we will concentrate on the large deviation regime, where one can 
expect 

with a positive AFj{qi,q2). The exponent L^''AFj{qi,q2) represents an interface free- 
energy cost to maintain the constrained values. In this regime, the large deviation 
functional can be expected to be self-averaging and one needs to compute the average 
free-energy 

E\ogiPj,Liqi,q2)) = -L^'AF(gi,g2). (4) 

The RSB implies that if one chooses qi = q2 in the domain where p{q) is non-zero, then 
AF{qi, qi) = 0. The property of homogeneity of the overlap should then translate in a 
form for AF(gi, ^2) of the kind AF(gi, ^2) = 9i{qi + '?2)/2) \qi - q2\^''- 

From the complete knowledge of the joint distribution ([3]) one can in principle 
extract the marginal distribution of the difference u = qi — q2 

Pl{u) = j dqi dq2 5 {u - {qi - ^2)) Piiqi, qi)- (5) 

In the large deviation regime, where this analysis is supposed to be valid, we have 
\qi ~ q2\ = 0{1) for L — > 00 and hence the integral should be dominated by the value 
of g = (gi + ^2)72 that maximizes the function g. As a consequence Pl{u), the PDF 
of u, should behave as e-3*^''l«l'^ in the tails. In order to understand the behavior 
for "small" values of u we can suppose a smooth cross-over between the small and the 
large fluctuation regimes. In this case the form (jll) suggests a flnite probability for 
|gi-g2| = ©(L"^"/^"), i.e., a scaling form Pi(gi, ^2) = H{{qi + q2)/2,{qi-q2)L^''^'''') and 
the marginal probability of the difference u = qi — q2, Pl{u) = P{uL^''/^'^)L^'^^^'' . This 
quantity Pl{u) is particularly easily accessible in numerical simulations and contains 
the the information on the ratio between the interface exponents. The comparison 
of this method with the study of the overlap correlation function [28], as well as the 
comparison of our analytic results in section 5 with pertubative calculations [27] will be 
not attempted in this paper. 
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3. Hierarchical Spin Glasses 

Here we introduce hierarchical spin-glass models on Dyson lattices. We will consider 
two families of models, a first one that is better suited for analytic studies and a second 
one that is more adapted to numerical simulations. In both families the spins are 
associated to the leaves of a binary tree, see fig. ([I]). The distance between two spins 
that, rising up in the hierarchy, meet after / branches is naturally defined as 2K The 
first family of models we would like to define, is the natural Spin-Glass generalization of 
the Ferromagnetic Dyson model ([9J). The Hamiltonian can be constructed iteratively, 
connecting the two non-interacting systems of 2^ spins {k G N) to form a composite 
system of 2'^"''^ spins in the following way: 

^k+li^l^ S2k+i] = Hf}[Si, 5*2*:] + i//^[S'2fc+i, 5'2fc+i] — 

^ 1,2'^ + ! 

i<j 

having defined the Hamiltonian for a system at the first level of the hierarchy 

H({S,,S,) = -^^JS^S,. (7) 

The couplings Jij are independent and identically distributed Gaussian random variables 
with zero mean and unit variance. The family is parametrized by the value of a that 
tunes the decrease of the strength of the interaction with the distance. The interaction 
strenght decrease as a power of the distance, the model is therefore a hierarchical 
counterpart of the one- dimensional spin glass with power-law interactions [T3] which 
has received attention recently [HI [151 [HI [HI [H]. As in this case, the model is defined 
for a > 1/2. Exactly at the value a = 1/2 it reduces to the infinite range model 
[19] . Notice that the sum of the squares of the interaction terms that couple the two 
subsystems ^t^rw E!<j 4 scales as 22('=+i)(i-'^). This is on the order of the volume 
for the mean- field value cr = 1 /2 and smaller than the volume for the nontrivial regime 
a > 1/2. 

The second family we consider is a dilute version of model (Q. We consider a 
Hamiltonian for 2^ spins with a fixed number of terms P = a2^[a > 1) 

p 

H^[Si, 5*2^] = — J^Sit^Sj^, (8) 

where the are taken as ±1 variables with equal probability, and we choose the 
interacting couples of sites [i^-,]^) independently term by term and in a way that, if 
the distance between i^^ and on the binary tree with 2^ branches is 2", we put the 
coupling with a probability given by 

P(z^,j^) =Afc 2"(^-2'^), (9) 

where Ak is a constant chosen such that it normalizes the probability. For a = 1/2 
the model reduces to the Viana-Bray spin-glass model [20j on Erdos-Renyi random 
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Figure 1. In our models the spins are associated to the leaves of a binary tree with k 
generation. In the first model all the spins interact with each other with a strength that 
decreases as a power of the distance. In the second case, each spin interacts directly 
with a small number of other spins at random, with a probability of interaction that 
behaves as a decreasing power of the distance. 

graph. We will refer to the two models defined above as model I ([6]) and model II (IE]), 
respectively. Both models should belong to the same universality class for the same value 
of 0", with respect to the critical and the low temperature properties. For a G [1/2, 1] the 
models have a finite temperature spin-glass transition. The critical point has a classical 
(i.e., Gaussian) character for a G [1/2,2/3], while it has a non-classical character for 
0" G [2/3, 1]. As observed in [181 Eij, diluted models as the one in are convenient in 
numerical simulations since the number of interactions for each spin does not grow with 
the system size. Spending an equal amount of computational effort we can therefore 
study much bigger sizes than for model I, hoping that the finite-size corrections are 
comparable in both models. 

4. Interfaces in hierarchical models 

Let us discuss some scenarios for the behavior of the PDF Pl{(1i,<12) corresponding to 
different possible physical situations. 

The simplest physical situation is the paramagnetic state where long-range order is 
absent. In this case we have a finite correlation length ^ and for large k one can expect 
qi and q2 to be sums of ~ L/,^ independent terms. The resulting probability distribution 
of qi and q2 is a product of two independent Gaussians with variance proportional to 

We can then consider a condensed phase with only two pure states. In this case, 
the space average value of the overlap can take two values ±g*. Small fiuctuations will 
still be Gaussian as in the paramagnetic case. Due to the symmetry of the model, where 
groups of spins at a given level of the hierarchy are on the same foot, one can expect that 
large fiuctuations where qi ~ q* and q2 ~ —q* (or vice versa) will imply a free-energy 
cost AF = A2('=+i)2(i— ) and have a probability e-/5^2('=+i)2(i"-)_ 

The most interesting possibility is a spin-glass phase with RSB. In this case one 
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has a zero mode in the free-energy associated to the existence of a couple of states with 
overlaps taking values in a finite interval [8j • If one chooses qi = q2 = Q one finds a broad 
distribution pk+i{q), which for large k is close to the limiting distribution p{q). Consider 
first a system in absence of interactions at the k + 1-th level. The two subsystems are 
independent and Pj^^iiqi, q2) = Pk{(li)Pk{(l2) ■ If the interaction is switched on, we can 
expect a free-energy cost AF(gi, ^2) = g{(li + q'2)2^^^~'^-*'^|q'i — q2\^''- The overlap interface 
exponent is just given by 6g = 2(1 — a) due to the fact that the total interaction strength 
sguarec? between the two parts scales as 2^^^^'^)^. The value of the exponent bg as well as 
the function g can be computed supposing RSB at the level k. The detailed calculation 
is quite involved and is presented in the next section; here we just give the net result, 
valid in the regime gi, ^2 ^ [—Qea, Iea]- Neglecting prefactors we could not compute, it 
reads: 

Pk+i{quq2) ~ e-''«''^+«^)/2)2'*^-^'"=^^'l''-'^^l^ (10) 

Notice the appearance of the k ^ 00 function p{q) in the exponent of (fTOj) . From this 
formula one can extract the conditional probability of the difference u = qi — q2 for a 
fixed value of the semi-sum g = (gi + g2)/2: 

P.+l(«|g)~e-^(«)2^'^""'''=^''l"l^ (11) 

Equation ( |TT|) summarizes our prediction for the hierarchical model. For a fixed value of 
the sum, the difference of the overlap is distributed according to the exponential of the 
cube, which is different from a naive Gaussian guess. The coefficient of the exponential 
is equal to the function p{q) which can be evaluated in independent measurements. In 
numerical simulations it might be complicated to collect sufficient statistics to condition 
u to the value of q. One can then turn to the unconditional distribution. In the large 
deviation regime, where u ~ 0(1), this should be dominated for large k and finite u by 
the value of q which maximizes (JTTi) . For functions p{q) as the one commonly met in 
spin-glass systems this is the value q = qEA which is the largest possible value of the 
overlap in the thermodynamic limit. 

The form ( fTTil suggests that the order of magnitude of the typical fluctuations is 
\u\ ~ 2^/'^(^~'^)('^+^). In this regime, the unconditional distribution of u involves the 
convolution of ( iTTl) with a presently unknown prefactor and cannot be computed. All 
the form ffTTl) tells us for the unconditioned distribution in this regime is that the finite 
volume distribution of u admits the scaling form 

Pk+i{u) = j dqP,^Mq)Pk+M = W^(22/3(i-'^)('=+i)|«|)22/3(i--)(^+i). (12) 

In order to compute the scaling function W the knowledge of the prefactor in (|TT1) would 
be necessary. We remark that the exponent in this scaling is a strong consequence of 
RSB theory, a naive guess would have suggested a Gaussian distribution with scaling 
variable 22(1— u^. 
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Qi + Qi 



(13) 



(14) 



5. Replica derivation of formula 

We analyze model 1 through the replica method. In order to compute the free-energy, it 
is natural to consider a recursion that relates the average partition function of n replicas 
Sf, a = 1, n, with fixed mutual overlaps Qab = ^ J2i=i ^i^i- Defining 

n \ l,n / 2*= 

^ a=l / a<b y i=l 

where Ej denotes the average over the disorder, we can write: 

Zfc[Q] = exp (^^22(1-)'= TrQ2^ J VQ,VQ2Zk-i[Qi]Zk-i[Q2]S (q ^ 

For integer n this is an exact relation based on the independence of the Hamiltonians of 
the sub-systems at level k. In principle, the thermodynamical properties of the system 
are encoded in this recursion and in its analytic continuations for n 0. For example it 
can be used to set up an epsilon expansion for the calculation of the critical indexes for 
cr e [2/3, 1] [22]. In this paper we just use f|T4|) to study the probability distribution f|TT]) . 
This can be done with the technique of constrained replicas, introduced and discussed 
at length in [7] . In order to consider constrained free-energies for two replicas with fixed 
overlaps one should fix some of the elements of the matrix Q^j, r = 1, 2 to the values 
of the constraints. Writing n = 2n' the constraint reads = Q^a+n' a — 1r fo^^ 

r = 1, 2 and a = 1, n'. Introducing the replica free-energy at level k for fixed qi and 
q2, Zi:[Q\qi,q2] = e~^^'"^*=['^l'?i''?2l we see that one needs in principle to compute: 

n' 

xU^iQia+n'-Ql)KQla+n'-Q2)- (15) 
a=l 

This form suggests that for large k, supposing the knowledge about -Ffc[Q], the integral 
over Qi and Q2 can be performed by the saddle-point approximation. Notice that the 
interaction term is sub-extensive and scales as 22(^^'^)('^+^), while the partial free-energies 
scale as the volume 2'^. The interaction term does therefore not contribute to the saddle 
point, and the maximization with respect to the matrices can be performed separately 
in each of the two sub-systems. 

The value of the interaction term at the saddle point determines the free-energy 
difference AF(gi,g2)- 

We now study the consequences of the hypothesis that there is RSB in the system. 
Specifically, we suppose that in absence of any constraints RSB is described by a 
continuous Parisi function qpix) taking values between the two extremes —qEA and qEA- 
In the constrained problem, correspondingly, each of the matrices is parametrized 
by two functions q'^{x) and p^(x) with x G [0, 1]. As implied by the analysis in [7] the 
free-energy in each of the two sub-systems is then independent of q^ , and the function 
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q^{x) and p^{x) can be directly related to the function qpix) of the unconstrained system 
by the relations 

qF{2x) X < Xi/2 ( fo \ ^ in 

I \ i / \ Qf{2x) X<Xr/2 

q[x) = < qr Xr/2 < X < Xr PxW = < , (16) 



I M y f I "i 

x) x> x\ y 

where Xr is the value of x such that qpixr) = qr- We then see that the free-energy 
difference from the unconstrained case is entirely due to the interaction term, which can 
be evaluated using the saddle-point value of the matrices Qi and Q2: 



Q1 + Q2Y /"^^ / gi(a^) + g2(x) y ^ fpi{x)+p2{x) 



.(17) 



Substitution of f|T6|) leads to 

rjn^ / Qi + Q2 j ^^f_2 f dxqp^xf + I {,q2 - qF{.x)){qF{x) - qi] 



'XI 

where without loss of generality we have supposed qi < q2- The first term is just the 
contribution that can be expected if qi = q2- Together with the two subsystems' free- 
energy it just gives the free-energy of the system at the level k + 1. The second term is 
associated to the free-energy excess needed to impose qi ^ q2- For small q2 — qi we can 
expand this last term and find 

nx' ( ^1 (g, - q,r, (19) 



2 

where x{q) is the inverse function of qrix). Inserted in (fTSll and identifying x'{q) with 
p{q) as discussed in length in [8] leads to f|TT]) 



6. Numerical Simulations 



In this section we discuss the results of numerical simulations used to test the behavior 
of the probability distribution of the overlap difference in the two subsystems. In order 
to deal with PDF's of a single variable we concentrated to the unconditional probability 
for which we have the theoretical prediction f|T2l) . We did not try to test the more 
detailed prediction ffTTl) which would need the numerical determination joint PDF's of 
two variables. 

We have simulated model II for a = 3/2, using parallel tempering to thermalize 
the system at low temperatures. 

We have tested our theoretical predictions concentrating on two values of a: 
cr = 0.576 which lies in the classical region 1/2 < a < 2/3, where the spin-glass 
transition is well described by mean- field theory, and a = 0.707 which lies in the non- 
classical region 2/3 < cr < 1, where the exponents are nontrivial. 
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Figure 2. Probability distribution function Pk{u) for K = 10, a — 0.576 (left) in 
the classical region and a — 0.707 (right) in the non-classical region, plotted in a log- 
linear scale as a function of u'^. Both systems are in the low temperature phase, the 
temperatures being respectively T — 0.745 and T = 0.740. The data for a = 0.576 
have been shifted of one unity to the right. For comparison we plot curves of the kind 
a e**" with a and b adjusted to fit the tails of the curves. 

We have characterized the critical point of the model following the procedure that 
we describe in the appendix. We estimate the critical temperatures to be equal to 
T, = 1.24(1) for a = 0.576 and = 1.14(1) for a = 0.707. 

To test the predictions described in the previous sections we measured the 
distribution of the variables q and function of temperature and system size. 

We considered systems sizes of 2'^, 2^, 2^ and 2^° spins. Averages were performed over 
4000 samples for the smaller systems and 700 samples for the largest system. The 
configurations are thermalized during the first 2^° Monte Carlo (MC) sweeps of the runs 
and then data is collected for the subsequent 2^° MC steps. 

The first prediction we test is the validity of the cube-exponential form in the tails 
for large u. This is well observed in all our simulations. A typical example of our 
findings is depicted in figure ([2]) where we plot the function Pk{u) for k = 10 for the 
two different values of a in the classical and non-classical regime, respectively. The data 
are plotted together with functions of form /(x) = a e"^'"^'^ which should be considered 
as a guide to the eye rather than the best fit. The parameters a and b were fixed by 
eye to be equal to a = 1, 6 = 40 for a = 0.707 and a = 0.8, 6 = 73 for a = 0.576. A 
best fit procedure results to be sensitive to the chosen fitting interval and to the tails 
of the distribution that represent probabilities too small to be correctly estimated with 
our statistics. Despite these caveats we believe that our data give an indication in favor 
of the cubic behavior for both values of a. 
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We then investigated the scahng of the PDF of u with the system size, contrasting 
it with the behavior of the PDF of q. In figures [3] and [5] we display the function Pk{q) 
for the two values of a at low temperatures. The function Pk{q) has the characteristic 
appearance of the one of systems developing RSB for large volumes, with two symmetric 
peaks and a non-zero part for q ^ 0. Figures H] and [6] in contrast show that the 
distribution of u, Pk{u), is unimodal around zero. Its width is as expected a decreasing 
function of k. In the insets the unsealed data are presented, while the in the main 
panel the result of scaling the data using the variable u2 shown. We judge the 

scahng (fT2l) to be very satisfactory, though not perfect for the values of k we explored. 
Indeed, while for a = 0.576 we could not find a the value of the ratio 9q/bq producing a 
better data collapse than the theoretical value |(1 — a) = 0.28, for a = 0.707 the value 
Oq/bq = 0.22 produces a better data collapse than the theoretical value |(1 — cr) = 0.19. 
We believe that this discrepancy is due to finite size effects, but we can not exclude at 
present that the theory should be amended in the non-classical region. 

Our data show that the best fitting exponent is largely temperature independent 
in the low temperature region. A direct scaling of the data in the high temperature 
region for the values of k = 7, 8, 9, 10 that we dispose, produces an effective exponent 
that crosses over slowly from the low temperature value towards the paramagnetic value 
Oq/bq = 0.5. We believe that we are seeing a preasymptotic behavior due to the influence 
of the critical fixed point. This influence could be particularly marked due to the 
power-law interactions where the critical fixed point continues to attract the system 
on relatively large scales. Larger and larger values of k are necessary to observe the 
paramagnetic behavior closer and closer to the critical temperature. 

7. Summary and Conclusion 

In this paper we have proposed a new method to determine numerically the overlap 
interface exponents first defined in [7]. These can be obtained by looking at the finite- 
size scaling of the PDF of the differences of the overlaps between two replicas in two 
subsystems. We have applied the definition to hierarchical models, where we could give 
a theoretical prediction for the scaling of the overlap differences with size. Though 
the interface exponent is naturally dictated by the model, the dependence in u is 
found to be nontrivial in presence of RSB. We tested this dependence in numerical 
simulations finding very satisfactory agreement both in the classical and in the non- 
classical region. These results confirm the interest of hierarchical spin-glass models, 
that combine analytical tractability, nontrivial critical points and RSB low temperature 
phases. 
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Figure 3. Probability distribution function p{q) for a = 0.576 (classical region) and 
T = 0.745 = 0.615Tc {T^ = 1.24(1)) for 1 = 2'' with fc = 7, 8, 9, 10. The peaks become 
sharper and sharper for larger k. 
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Figure 4. Scaling plot of probability distribution of the overlap difference for 
a = 0.576 (classical region) and T = 0.745 = 0.601Tc for L = 2^ with k = 7,8,9, 10. 
Inset: the unsealed probability. 
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Figure 6. Same as figure|3]for a = 0.707 in the non-classical region. The temperature 
is the same as in figure [3 T = 0.740. 
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Appendix: Characterization of the critical point. 

In this appendix we discuss the characterization of the critical point of the model for 
the values of the interaction parameter a that we have considered in the text. 

The model has been simulated with the Parallel Tempering algorithm, with 10 
values of the temperature, using 2^° thermalization steps before collecting data in the 2^° 
steps. We considered systems sizes of 2'', 2*, 2^ and 2^° spins. Averages were performed 
over 4000 samples for the smaller systems and 700 samples for the largest system. 
Simulating two replicas in parallel, we have studied the second and the fourth moment 
of the distribution of the mutual overlap q2 = E{Q^) and g4 = E{Q^). We have identified 
the critical temperature and the exponents t] and u using finite size scaling through the 
behavior of q2, and the corresponding Binder parameter [23] S = |(3 — ^). In the 
non-classical region, cr > 2/3 where finite-size scaling should hold the various parameters 
exhibit the following dependence on the temperature and system size L = 2^: 

X = Lq2 = L''^g2{L'/''{T-T,)) (21) 
B = b{L^/''{T ~T^)). 

The exponent rj should not renormalize in long range models, and analogously to the 
Euclidian ID model take the value r] = 3 — 2a both in the non-classical and in the 
classical regions [13j. 

In the classical region, 1/2 < a < 2/3, the scaling implied by ( 1211) does not hold 
[211 125]. It is possible to show, using the fact that the critical theory is described by 
a cubic action analogous to the one for short-range spin glasses [26], that the various 
quantities scale according to the following: 

Lq2 = \T- T,\-^~g2{L^'\T - T,)) (22) 
B = b{L^/^{T -T,)). 

The exponent 1/3 in the scaling functions can be derived from dimensional analysis 
from the cubic action. The exponent 7 takes the value jci = 1 independently of a as 
can be checked from the scaling relation 7 = (2 — ri)v with v = 2^zi- 

Let us now turn to the data considering the non-classical region first. To analyze the 
data we observe the following procedure: we first estimate the critical temperature from 
the crossing point of the Binder parameter and the rescaled values, especially, ^2-^^^^ 
and g4L~^+^''' since these provide for a cleaner crossing than the Binder parameter. We 
then fix the value of u in order to collapse the curves. The result for cr = 0.707 is 
shown in figure 13, we present the data for q2, q^, and B scaled as in fl22]) using a value 
of = 1.14. 

If we try to use the scaling (121 p in the classical regime we get inconsistent results: 
although we obtain an approximate crossing of the curves for the three quantities, the 
temperatures at which the curves cross do clearly not coincide. The crossing of the 
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Figure 7. Behavior of q2, and _B as a function of temperature and system size 
for cr = 0.707. In tlie inset we plot 92/-^^"'', 94/-^^"^'' and B against (T ~ TJ, with 
the values Tc — 1.14 and 77 = 3 — 2a. Main panel, scaling plot of the data q2/L^~'^, 
qi/L^-'^'^ and B against (T - T^)L^/'' assuming the value oi v = 3.09. 



Binder parameter indicates T^. = 1.24 and as shown in figure [7| for a = 0.576, we obtain 
consistent scaling assuming tlie form fl22l) and not fl2T]) . 
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Figure 8. Behavior of q2, qi and _B as a function of temperature and system size 
for (T = 0.576. In tlie lower inset we show the Binder, for which the data collapse 
is obtained for Tc = 1.14. The main panel shows Lq2\T — Tc\- Notice the high 
temperature behavior corresponding to the Curie law. The curves for larges systems 
sizes show a tendency of converging toward the value 1 close the critical temperature. In 
the upper inset we show the scaling for L^(74|T — which shows the same qualitative 
behavior. 
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